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ABSTRACT 

The nucleosynthesis of Be and B by spallation processes provides unique insight 
into the origin of cosmic rays. Namely, different spallation schemes predict sharply 
different trends for the growth of LiBeB abundances with respect to oxygen. "Primary" 
mechanisms predict BeB oc O, and are well motivated by the data if 0/Fe is constant 
at low metallicity. In contrast, "secondary" mechanisms predict BeB oc and are 
consistent with the data if 0/Fe increases towards low metallicity as some recent 
data suggest. Clearly, any primary mechanism, if operative, will dominate early in 
the history of the Galaxy. In this paper, we fit the BeB data to a two-component 
scheme which includes both primary and secondary trends. In this way, the data 
can be used to probe the period in which primary mechanisms are effective. We 
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analyze the data using consistent stellar atmospheric parameters based on Balmer 
line data and the continuum infrared flux. Results depend sensitively on Pop II O 
abundances and, unfortunately, on the choice of stellar parameters. When using 
recent results which show 0/Fe increasing toward lower metallicity, a two-component 
Be-0 fits indicates that primary and secondary components contribute equally at 
[0/H]gq = -1.8 for Balmer line data; and = -1.4 to -1.8 for IRFM. We apply 

these constraints to recent models for LiBeB origin. The Balmer line data does not 
show any evidence for primary production. On the other hand, the IRFM data does 
indicate a preference for a two-component model, such as a combination of standard 
GCR and metal-enriched particles accelerated in superbubbles. These conclusions rely 
on a detailed understanding of the abundance data including systematic effects which 
may alter the derived 0-Fe and BeB-Fe relations. 
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1. Introduction 

The light elements beryllium and boron (BeB) provide unique insight into the nature 
of non-thermal nucleosynthesis in our Galaxy. Due to their low binding energy, these nuclei 



are not produced significantly in the big bang (Thomas, Schramm, Olive, & Fields 1993 , 



Delbourgo-Salvador & Vangioni- Flam |199|) or in stellar nuclear burning. Instead, BeB are made 



by spallation processes due to energetic nuclei and neutrinos. The nucleosynthesis origin of BeB is 
in principle encoded in their abundances in the most primitive, metal poor (Population II) stars. 
Indeed, the mere presence of Be in halo stars is perhaps the strongest evidence that accelerated 
particles were present in the early Galaxy.]^ Thus, in this paper we will focus on Pop II BeB 
trends, via a detailed analysis of the data and a comparison with recent models. 

The production of BeB by Galactic cosmic ray (GCR) spallation of interstellar CNO nuclei 
was the standard model for BeB nucleosynthesis for almost two decades after first being proposed 
(Reeves, Fowler, & Hoyle [1970 ; Meneguzzi, Audouze, & Reeves 1971| ). However, this simple model 



was challenged by the observations of BeB abundances in Pop II stars, and particularly the BeB 
trends versus metallicity. Measurements showed that both Be and B vary roughly linearly with 
Fe, a so-called "primary" scaling. In contrast, standard GCR nucleosynthesis predicts that BeB 
should be "secondary" versus spallation targets CNO, giving Be oc O^ (Vangioni-Flam, Casse, 
Audouze, & Oberto [T990| ). If O and Fe are co-produced (i.e., if 0/Fe is constant) then the data 



clearly contradicts the canonical theory, i.e. BeB production via standard GCR's. 

BeB abundances in Pop II stars have thus provided the motivation for primary BeB 
nucleosynthesis, in which the production rate is independent of metallicity. These primary models 
invoke the extraction and acceleration of fresh products of nucleosynthesis (a's, CO) originating 
from Type II supernovae and Wolf-Rayet stars and fragmenting on interstellar H and He. They 
should be considered in addition to the standard GCR (secondary) process; any primary process 
would be expected to play a major role in the early Galaxy (halo phase), whereas the secondary 
one takes over in the galactic disk (Casse et al. 1995, Vangioni-Flam et al. 1996) or earlier as we 
shall see. The transition between the two modes may vary depending on the adopted BeB data, 
as well as the [0/Fe] vs [Fe/H] correlation. 

One such primary model has focused on superbubbles-regions of hot, rarefied, metal-rich gas 
swept out by the collective effects of massive star winds and supernova explosions. These regions 
have been proposed (Vangioni-Flam et al. 199^ ; Higdon et al. 1998| ; Ramaty & Lingenfelter 1999; 



Parizot 1998; Vangioni-Flam, Casse & Audouze 1999) to accelerate freshly synthesized, metal-rich 
particles. The energy spectrum of these nuclei has been extensively studied by Bykov (1999 jl995| ) 
and it is remarkably similar to the GCR injection spectrum. Finally, core collapse supernovae 



^Pop II lithium abundances are dominated by the primordial ^Li component, which must be subtracted to provide 



information about Galactic Li sources. Elemental Li is discussed in Ryan, Beers, Olive, Fields, & Norris (L999); the 
Li isotopes are analyzed in Vangioni-Flam et al. ( |L999| ) and Fields & Olive ( |L999b| ). 
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themselves have been proposed as a site of ^^B production (and some via spahation reactions 
between supernova neutrinos passing through the ^^C layer, the 'V-process" (Woosley et al. 1990). 



Recently, another solution has been proposed to resolve discrepancy between the observed 
BeB abundances as a function of metallicity and the predicted secondary trend of GCR spallation. 
As noted above, standard GCR nucleosynthesis predicts Be oc O^, while observations show 
Be ~ Fe, roughly; these two trends are inconsistent if 0/Fe is constant in Pop II. However, recent 



observations find 0/Fe increasing at low metallicities (Israelian, Garci'a-Lopez, & Rebolo 1995; 



Boesgaard et al. |1999a|) . In fact Boesgaard et al. ( |1999aD argue for a single slope of ~ —0.35 



for 0/Fe vs. Fe at all metallicities. As shown by Fields & Olive ( |l999a ), a combination of 



standard GCR nucleosynthesis, and z^-process production of ^^B is consistent with current data. 
One should note that while the trend in 0/Fe is seen in the OH data and agrees with the O I 
triplet (Boesgaard et al. |1999a ), they continue to disagree with the observations of 0/H using the 



forbidden [O I] hue (Fulbright & Kraft 1999). 

It is clear that spallation processes, BeB abundances, and 0-Fe evolution are closely linked. 
In this paper, we explore these connections via careful analysis of both data and theory. We 
use the most current BeB data, with abundances derived from consistent stellar atmosphere 
models. We consider stellar parameters based the Balmer-line data of Axer, Fuhrmann, & Gehren 
(|1994f ) and the IRFM (Alonso et al. |1996bD . We first fit the data to both iron and oxygen to 



determine the overall primary vs. secondary nature of the BeB abundances. These data are fit to 
a two- component, primary plus secondary, metallicity dependence. This fit quantifies both the 
strength of each component as well as the metallicity and hence epoch as which the secondary 
component becomes dominant. 

We compare these phenomenological results to the predictions of current models for BeB 
nucleosynthesis and chemical evolution. We find that for the Balmer line data, any primary 
component could only have been dominant at oxygen abundances [0/H] < —1.8, which is at 
metallicities below the existing data. That is, for this data set, we find no evidence for primary BeB 
production. In contrast, for the IRFM data, we find that primary production dominates for [0/H] 
< —1.4 to —1.8 (depending on the exact data-see below). The higher value indicates indicating 
the need for a primary mechanism at low metallicities. If 0/Fe is constant then the transition 
between primary and secondary production processes would occur at [Fe/H] ~ — 1. As such, our 
quantitative conclusions can not be definitive. The derivation of the oxygen abundance from the 
observed spectra is delicate, specifically at low metallicity. Moreover, non local thermodynamic 
equilibrium (NLTE) effects on Fe could affect substantially the BeB-Fe relationship. We also 
point out that since magnesium is not plagued by the same difficulties regarding the stellar yields 
of Fe, it would be worthwhile to examine the relationship of the BeB elements with Mg. New 
observations are eagerly awaited. 
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2. Data Sets 

The BeB-OFe abundance trends encode the history of spallation in the Galaxy. Our ability 
to infer this history is completely determined by the accuracy of the abundance data studied. 
It is thus necessary to obtain high quality data. Just as importantly, the abundances must be 
combined in a systematic and consistent way. 

The determination of abundances from raw stellar spectra requires stellar atmosphere models. 
The atmospheric models require key input parameters, notably the effective temperature Tgfj 
and surface gravity g, and assumptions regarding, e.g., the applicability of local thermodynamic 
equilibrium (LTE). Unfortunately, there is no standard set of stellar parameters for the halo 
stars of interest. In practice, different groups derive abundances via different procedures, which 
give similar results but retain systematic differences. The systematic differences in the data can 
in fact obscure the BeB-OFe trends one seeks. Thus, to derive meaningful BeB fits, one must 
systematically and consistently present abundances derived under the same assumptions and 
parameters for stellar atmospheres. 

Below, we will present results based for the available BeBOFe data based on two methods 
of analysis. We will refer to these as the Balmer line data and the IRFM data. The Balmer 
line data is taken from the work of Axer, Fuhrmann, & Gehren ( |1994 ) and Fuhrmann, Axer, 
& Gehren ( |1993|) . These authors have studied over 100 dwarf and subgiant stars to determine 
a consistent set of stellar parameters which includes the effective temperature, surface gravity 
and iron abundance. The effective temperatures are determined from the synthesis of the first 
four Balmer lines. This amounts to a spectroscopic rather than a photometric determination of 
the effective temperature. The remaining parameters were determined simultaneously from the 
analysis of a number of iron lines. The temperatures determined this way are typically 100-200 
K higher than many of the photometric determinations in the literature. The [Fe/H] abundance 
determined by Axer, Fuhrmann, & Gehren ( |1994| ) are also typically higher than others sometimes 
by as much as 0.5 dex. It goes without saying that such differences can have an enormous effect 
in trying to establish a trend for BeB vs. Fe. 

An alternative method to obtain effective temperatures is achieved by comparing the 
bolometric flux and the IR flux at a fixed wavelength. This method, knows as the IRFM, 
requires only the theoretical prediction of the continuum IR flux. There is a slight dependence 
on the assumed surface gravity and metallicity. A large compilation of nearly 500 stars for 
which the IRFM was used to determine temperatures was presented in Alonso, Arribas, & 
Martinez- Roger ( |1996a| ). Bonifacio and Molaro used these temperatures ( [L997| ) to accurately 
determine Li abundances in 41 plateau stars as a follow-up of their previous work using Balmer 
line temperatures (Molaro, Primas, & Bonifacio 1995| ). For Li similar results were found for the 



plateau abundance in each of the two methods. In Alonso, Arribas, & Martinez- Roger ( |l996b ). 



calibrated expressions for the effective temperatures were derived based on the earlier tabulation. 
Unfortunately, these calibrations do not always lead to the same temperatures found in the 
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tabulations and differences by as much as 200 K occur. We will present results based on both (1) 
the IRFM temperatures in the compilation Alonso, Arribas, & Martinez- Roger ( p.996a ), and (2) 
based on the calibration. 

We can not overly stress the importance of reliable stellar data. The Balmer-line data appear 
to be self-consistent, and are probably the most reliable. However, because we will present results 
based on the IRFM temperature scales, we would like to point out that there are significant 
differences in the data. To illustrate the point we take for example the case of the star BD 3° 
740. From Axer et al. (|1994|) , we find this star to have {T^sMg, [Fe/H]) = (6264, 3.72, -2.36). 
The beryllium and oxygen abundances for this star was reported by Boesgaard et al. ( |l999b| ) 
and ( 1999a ). When adjusted for these stellar parameters, we find [Be/H] = -13.36, and [0/H] 
= -1.74. In contrast, the stellar parameters from Alonso et al. (|1996a| ) are (6110,3.73,-2.01)0 
with corresponding Be and O abundances of -13.44 and -2.05. Garcia-Lopez et al. ( |199^ ) use a 
calibrated IRFM based on Alonso et al. ( [1996b| ) and take (6295,4.00, -3.00). For these choices, we 
have [Be/H] = -13.24 and [0/H] = -1.90. Notice the extremely large range in assumed metallicities 
and the difference in the two so-called IRFM temperatures. While this may not be a typical 
example of the difference in stellar parameters, it is differences such as this (and this star is not 
unique) that accounts for the difference in our results and the implications we must draw from 
them. 

We also note that in addition to the large systematic differences in the Fe abundances (an 
order of magnitude in the above example), it has recently been argued that NLTE effects for Fe, 
particularly in low metallicity stars are non-neglig ible (Thevenin & Idiart |1999| ). At the low end 
of the metallicity range considered here, they argue for a ~ 0.3 dex upward correction in [Fe/H]. 
We have not included this correction in the data discussed below. 

In the next section we will fit the BeB-OFe abundances over Pop II metallicities. The fits span 
a metallicity range —3 < [Fe/H] < —0.8, and —2 < [0/H] < —0.5 depending on the particular 
choice of stellar parameters. The upper bounds roughly mark the disk-halo transition, and the 
low metallicity bounds are just set by the availability of BeB data. 

As we indicated, we will consider three sets of stellar parameters to be used in the BeB vs OFe 



analysis. The three are: stellar parameters based on the Balmer line data of Axer et al. (|1994| ) 
to be denoted as Balmer; the IRFM, with data from Alonso et al. ( 1996a| ) - denoted as IRFMl 



the IRFM with stellar parameters as determined from the calibrations in Alonso et al. (1996b) 
as reported by Garcia-Lopez et al. ( [1998| ) and Israelian et al. ( 1998 ) - denoted as IRFM2. Note 
that at low metallicities this calibration is based on an analytic formula which is divergent. Thus, 
the points at low metallicity could be questionable. Be abundances used in the present work are 
from Rebolo et al. ([T988D , Ryan et al. ( [T990i p92D , Gilmore et al. (|l99^ ), Rebolo et al. ( p9^) , 
Boesgaard & King ( pJ9^ ), Primas (|l99|), Hobbs & Thorburn (|l99|), Molaro et al. ( |l997| ). 



^The surface gravity was recalculated by Bonifacio & Molaro (1997) 
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Boesgaard et al. ( |1999b[ ). For each star, and for each choice of Balmer, IRFM1,2, we first adjust 
the Be abundance to a common set of stellar parameters. Then when multiple observations for 
a given star are available these abundances are weight averaged. LTE is commonly used, since 
non-LTE (NLTE) corrections for Be are not expected to be significant. In total, we have data on 
Be for 37 stars. Not all of the stellar parameters have been determined in a uniform way for the 
entire set of 37 stars. In the case of the Balmer set, we have parameters from Axer et al. ( |1994| ) for 
22 stars. This set can be enlarged somewhat using the data of Boesgaard et al. ( |l999bD who have 
presented results based on two scales, the King ( |1993| ) scale which they claim is similar to Balmer 
and Carney ( 1983| ) which they claim is similar to IRFM. Thus for Balmer + King we have data 
for 30 stars. In the case of IRFMl, we also have data on 22 stars. IRFM 1 + Carney contains 25 
stars. Finally for IRFM2, there are only 18 stars, but IRFM2 + Carney gets us up to 24 stars. 

B abundances are taken from Duncan et al. ( |1997| ) and Garcia-Lopez et al. (1998) and Primas 
et al. (1999). NLTE corrections for B are significant (Kiselman |1994 , Kiselman & Carlsson 1996), 
and are not uniform over metallicity, giving larger enhancements at the lower metallicities. Thus, 
NLTE correction have the systematic effect of flattening the B-OFe trends. We apply the same 
procedure as described above for the boron data. There are a total of 15 low metallicity stars with 
boron measured, of which we have Balmer data for 11 (and Balmer + King on 13). 9 stars are 
known for IRFMl (and 10 for IRFMl + Carney), and 11 for IRFM2. 



Table 1: Slopes for 0/Fe versus Fe. 



number 


method 


slope 


22 


Balmer 


-0.48 ±0.16 


31 


Balmer 


-0.51 ±0.10 


36 


Balmer 


-0.47 ±0.09 


22 


IRFMl 


-0.45 ±0.14 


27 


IRFMl 


-0.44 ±0.10 


36 


IRFMl 


-0.43 ± 0.09 


22 


IRFM2 


-0.35 ±0.10 


29 


IRFM2 


-0.36 ±0.09 


36 


IRFM2 


-0.32 ± 0.08 



Pop II oxygen abundances, and 0/Fe, are critical for all studies of nucleosynthesis. The 
long discussion of O abundances in the literature illustrates that the path from stellar spectra 
to abundances is not a trivial one. Three different lines have been used to obtain the oxygen 
abundances with varying results. These are: (1) the allowed O I triplet at 7774 A; (2) the forbidden 
[O I] line at 6300 A; and (3) molecular OH lines at 3085 A. The third method was suggested by 
Bessell, Hughes & Cottrell (|1984 ) and recently high resolution UV spectra were obtained for many 
of the stars with Be and B observations (Nissen et al. 1994, Israelian, Garcia-Lopez, & Rebolo 
|1998| , Boesgaard et al. 1999a). These new data suggest a significant nonzero [0/Fe] trend (versus 
[Fe/H]) in Pop II stars. We use all of the recent data on OH lines and combine the data in a 
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similar manner as that described for Be. In the two most recent works, a strong case is made for 
a varying 0/Fe at low metallicity. The results of these two groups are in excellent agreement. It 
should be noted that the forbidden line continues to show lower O/H abundances as claimed in 
the recent work of Fulbright & Kraft ( |199S| ) . 

For O/H, we have a total of 36 stars with low metallicity OH data. There are 22 stars with 
Balmer data, and 31 with Balmer + King. There are 22 stars with IRFMl data and 27 with 
IRFMl + Carney. There are also 22 IRFM2 O/H data points and 29 IRFM2 + Carney. We wih 
define the logarithmic slope WQ/Pe ^^r 0/Fe vs Fe by 

[0/Fe] = a;o/Fe[Fe/H] + const (1) 

Our results for the slopes are given in Table ||. Within the uncertainties these slopes are all 
consistent with one another. The IRFM2 slopes are somewhat smaller and in very good agreement 
with those found by Israelian et al. ( 1998| ) as should be expected since the IRFM2 temperature 
scales are taken from that work. 



3. Statistical Tests of BeB Origin 

We will first derive the overall fits for BeB vs both O and Fe using the data discussed above. 
The Be and B data can be fit versus [Fe/H] or [O/H], using logarithmic abundances so that for 
example, 

[Be] = WBeFe [Fe/H] + const (2) 

We will focus on the logarithmic slopes such as t^BeFe; the same procedure, applied to boron, gives 
LVBFe, etc. [Be] is defined as log(Be/H) + 12. Our results are summarized in tables below. 

Note that there is a sharp difference in slopes of BeB vs Fe/H as compared with O/H. This is 
implied by the variability in 0/Fe at low metallicity. Also note that for the cases of Balmer and 
IRFMl, the scaling of Be with respect to O/H is almost purely secondary. Boron, on the other 
hand shows a strong primary component, which in this case, we would expect to be due to the 
z^-process (Woosley et al. 1990, Olive et al. 1994| ). The IRFM2 data still indicates the need for an 



early primary Be component as we explain below. 

We would also like to call attention to the slopes for B/Be with respect to either Fe/H or O/H. 
They are all non-zero (even when the relatively large uncertainties are taken into account), though 
there is a considerable dispersion in the data. This is a departure from the standard picture of 
purely secondary GCR nucleosynthesis as well more recent primary models of BeB nucleosynthesis. 
This behavior is due in part to the new Be observations of Boesgaard et al. ( |l999bD for some 
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Table 2: Slopes for Be versus Fe and O. 



method 


number 


tracer 


slope 


number 


tracer 


slope 


Balmer 


22 


Fe 


1.39 ±0.16 


19 





1.78 ±0.19 


Balmer 


30 


Fe 


1.26 ±0.11 


27 





1.79 ±0.16 


Balmer 


37 


Fe 


1.27 ±0.10 


31 





1.70 ±0.17 


IRFMl 


22 


Fe 


1.23 ±0.14 


21 





1.83 ±0.19 


IRFMl 


25 


Fe 


1.19 ±0.11 


24 





1.80 ±0.17 


IRFMl 


37 


Fe 


1.21 ±0.10 


31 





1.49 ±0.14 


IRFM2 


18 


Fe 


1.18 ±0.11 


18 





1.36 ± 0.09 


IRFM2 


24 


Fe 


1.15 ±0.10 


24 





1.35 ±0.08 


IRFM2 


37 


Fe 


1.18 ±0.09 


31 





1.29 ±0.08 



Table 3: Slopes for B versus Fe and O. 



method 


number 


tracer 


slope 


number 


tracer 


slope 


Balmer 


11 


Fe 


0.78 ±0.22 


10 





1.23 ±0.32 


Balmer 


13 


Fe 


0.56 ±0.14 


12 





1.22 ±0.28 


Balmer 


15 


Fe 


0.60 ±0.14 








IRFMl 


9 


Fe 


0.73 ±0.19 


9 





0.98 ± 0.28 


IRFMl 


10 


Fe 


0.63 ±0.13 


10 





0.94 ±0.23 


IRFMl 


15 


Fe 


0.70 ±0.10 


12 





1.00 ±0.20 


IRFM2 


11 


Fe 


0.72 ±0.14 


11 





1.02 ±0.16 


IRFM2 


15 


Fe 


0.70 ±0.11 


12 





1.00 ±0.16 



relatively low metalicity stars. There are now a few cases for which the B/Be ratio in fact exceed 
100 as can be seen in the figures. This is true for all of the choices of stellar parameters considered. 



This in fact is a prediction of the secondary model with z^-process B production (Olive et al. 1994 
Fields, Olive, & Schramm 1995). 



We are now in position to test the BeB data for the presence of two components, primary 
and secondary. Because physically, BeB production is tied to elements such as CNO, for the data 
available, we analyze only the Be and B trends versus oxygen. For A G BeB, the data are fit to 



A 
1 



O 

Ol 77 ± 02 



H 
A 
H 







O^ 
0/H 

(oTh)^ 



/ O/H 



(3) 
(4) 



where the fit parameters Oj are expressed in "scaled" units ai = {0/A)Qai and 
02 = (0/11)0/(74/11)002. Note that in eq. (|^, the abundances are not logarithmic, but 
are the (linear) ratios with respect to hydrogen. The fit parameters ci and 02 (equivalently, ai 
and 02) quantify respectively the primary and secondary contributions to A. These are given in 
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Table 4: Slopes for B/Be versus Fe and O. 



method 


number 


tracer 


slope 


number 


tracer 


slope 


Balmer 


9 


Fe 


-0.52 ±0.24 


9 





-0.68 ±0.37 


Balmer 


11 


Fe 


-0.57 ±0.16 


11 





-0.91 ±0.30 


IRFMl 


9 


Fe 


-0.49 ±0.18 


9 





-0.81 ±0.27 


IRFMl 


10 


Fe 


-0.59 ±0.15 


10 





-0.96 ±0.25 


IRFMl 


11 


Fe 


-0.60 ±0.15 


11 





-0.71 ±0.20 


IRFM2 


10 


Fe 


-0.57 ±0.14 


10 





-0.53 ±0.21 


IRFM2 


11 


Fe 


-0.56 ±0.14 


11 





-0.49 ± 0.20 



Table |. 



Table 5: Break points for Be versus O. 



number 


method 


ai 


a2 


[0/H]eq 


19 


Balmer 


0.042 ± 0.003 


2.30 ± 0.70 


-1.75 


27 


Balmer 


0.027 ±0.027 


2.21 ±0.49 


-1.94 


31 


Balmer 


0.027 ±0.027 


2.28 ±0.52 


-1.94 


21 


IRFMl 


0.034 ± 0.034 


2.14 ±0.53 


-1.79 


24 


IRFMl 


0.031 ± 0.027 


2.28 ±0.51 


-1.88 


31 


IRFMl 


0.050 ± 0.031 


2.14 ±0.54 


-1.62 


18 


IRFM2 


0.111 ±0.031 


2.57 ±0.76 


-1.37 


24 


IRFM2 


0.115 ±0.031 


2.61 ±0.68 


-1.36 


31 


IRFM2 


0.122 ±0.031 


2.79 ±0.70 


-1.35 



The terms in eq. (^) are equal at the metallicity 

with the primary term dominating at 0/H < (0/H)eq, and the secondary term dominating at 
0/H > (0/H)eq. Physically, (0/H)eq identifies the epoch at which primary cosmic ray sources 
are overcome by secondary sources; we thus refer to this as the "break point" of the A-0 trend. 
These points (in [0/H]) appear as [0/H]eq in Table ^, and indicate the metallicity at which the 
transition from primary to secondary should occur. Recall that any primary component present 
will be dominant at sufficiently low [0/H]. Therefore we can use the information in the table to 
compare the relevant strengths of the primary and secondary components and constrain models 
for cosmic ray origin and acceleration. As one can see, for both Balmer and IRFMl, the break 
point occurs at very low [0/H], in fact at the edge of data as can be seen from the figures where 
the data is plotted. 

There is also the question regarding how meaningful these two-parameter fits are. To give a 
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quantitative answer to this question requires an examination of the \^ of the various fits in Table 
^. In Table ^ we show the relevant x'^'s for the break point fits together with those for fits where 
either (X\ or a2 are forced to be 0. That is the for either a purely linear or purely quadratic 
evolution. 



Table 6: foi' the break points in Table 5. 
number method x^ X^(q^i = 0) (Pp) = 0) (Pp) 



19 


Balmer 


18.4 


20.2 (78%) 


30.9 (99.e 


%) 


27 


Balmer 


27.2 


28.3 (68%) 


54.1 (99.9 


9%) 


31 


Balmer 


34.1 


35.3 (68%) 


61.2 (99.9 


9%) 


21 


IRFMl 


13.8 


15.0 (68%) 


30.7 (99.9 


9%) 


24 


IRFMl 


18.1 


19.4 (78%) 


39.6 (99.9 


9%) 


31 


IRFMl 


34.4 


39.1 (94%) 


53.1 (99.9 


S%) 


18 


IRFM2 


24.2 


56.4 (99.97%) 


46.6 (99.8 


%) 


24 


IRFM2 


29.5 


61.8 (99.99%) 


54.0 (99.9 


7%) 


31 


IRFM2 


40.0 


72.7 (99.99%) 


65.7 (99.9 


3%) 



In order to justify an additional parameter, we must find a significant drop in the x^- The 
degree of significance can be quantified by the F-test, which produces a likelihood that the 
2-parameter fit is a better description of the data than the 1-parameter (primary or secondary) 
fit. This quantity is also shown in the Table ^ as Pp- For example, in the first line of that table, 
for the Balmer case with 19 points, we see that while the 2-parameter fit is clearly superior (with 
a probability of 99.6%) to a straight primary (linear) fit, it is only an improvement at the 78% 
level over a purely secondary fit. That is, in the latter case, it should be questioned whether or 
not the 2-parameter fit is actually better than a purely secondary (quadratic) fit. This pattern is 
similar for the IRFMl data set as well. In contrast, for the IRFM2 data, the 2-parameter fit is 
clearly preferred to either the purely primary or secondary fits. 

Before concluding this section, we note that use of Fe abundances to study the BeB evolution 
presents us with several difficulties. In addition to the overall uncertainties in the observational 
determination of [Fe/H] (with or without NLTE effects), we will also find it difficult to model 
the recently observed 0/Fe slopes. In section 4.4, we note that most calculated yields predict 
an overabundance of Fe with respect to O, particularly for large stellar masses, when compared 
with the recent observations. There is however, a substantial uncertainty in the ejected yield of Fe 
relative to the mass of Fe trapped in the stellar remnant. For this reason, it has been suggested 
(Shigeyama & Tsujimoto 199^ ) that the use of Mg instead of Fe should be a reliable tracer 



for chemical evolution studies. Indeed, there have been some recent preliminary observational 
attempts (Fuhrmann et al. 1997, Fuhrmann 1998) to obtain Mg data for the halo stars of interest. 



In this context, it happens that Mg is less affected by spectroscopic uncertainties than O; as such 
Mg might be a better metallicity tracer than either O and Fe. 
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4. Models 



Using the results of the BeB fits, we can now test different scenarios for their non-thermal, 
spallative nucleosynthesis. To compute the BeB-OFe trend expected for each scenario, one must 
combine models for cosmic ray /accelerated particle origin, particle propagation, and chemical 
evolution. We describe these in turn. 



4.1. Cosmic Ray Sources 



We focus on two suggestions for the origin of accelerated particles: 

1. Standard GCR, which leads to secondary Be and B trends (proportional to O^). 

2. Superbubble accelerated particles (SAP). The metal-enriched composition of the superbubbles 
is reflected in the particles and leads to primary Be and B production (proportional to O). 



For each of these, we must specify the source spectrum Q{E) = dN/dE of the accelerated particles, 
their composition, and the dependence of the composition on that of the ambient medium. For 
further discussion see Vangioni-Flam, Casse, & Audouze (1999). 



4.1.1. Standard Galactic Cosmic Rays 

For many years it has been a common perception that supernova remnants are the principal, 
perhaps predominant, agents which accelerate galactic cosmic rays up to an energy of around 10^^ 
eV. Supernova remnants have the energy necessary to satisfy the cosmic ray energetics requirement 
if their acceleration efficiency is about 0.1; the current estimate of supernova rates in our Galaxy 
can adequately supply the observed cosmic ray density. Such cosmic rays are presumed to be 
energized by diffuse shock (Fermi) acceleration, the shock itself being induced by the impact of 
the supernova ejecta on the surrounding interstellar medium (e.g., Ellison, Drury, & Meyer 1997| ). 



We follow here the standard approach, which assumes that supernova blast waves create 
strong shocks which accelerate ISM material, with a composition that retains certain biases. 
Specifically, the model follows very closely that of Fields & Olive ( [1999a ), but is also similar to 



that of Lemoine et al. ( |1998D . Since the acceleration engines are supernovae, the source intensity 
scales as the supernova rate and thus (for all but the earlies times) the star formation rate: 
Q{t) oc Nii{t) ~ tpit). Note that in the early Galaxy, Type II supernovae dominate. 

The source spectrum is that of particles accelerated by strong shocks. To a good 
approximation (e.g., Ellison, Drury, & Meyer |1997D , this gives a power law in momentum, of the 
form dN/dp oc p^^. In energy space, we have Q{E) = dN/dp dp/dE, which goes to q{E) ~ E^^'^ 
at low energies, and q{E) ~ E~'^ at high energies. Of course, the full, relativistic p{E) expression 
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is used for all energies in the numerical work. The proton source spectrum appears in Figure [T|a. 
This propagated spectrum (see Figure lb) is consistent with the one determined on the basis of a 
refined propagation model fitting a large number of observations by Strong & Moskalenko ( 19981 ). 



The accelerated particle composition in species i reflects that of the ISM: 

yfl^(t) = Q,{t)/Qp{t) a X,{t) 



(6) 



with the following scaling. The present-day observed cosmic-ray source composition, (to), 
differs from a solar (and ISM) composition in that refractory and/or low first ionization potential 
elements are enhanced with respect to the volatile and/or high ionization potential ones (Meyer, 
Ellison, & Drury 1997, Casse and Goret 1978, Vangioni-Flam and Casse 2000| ). For example, 
oxygen has a cosmic ray abundance y^'^ = 3.5 x 10"'^ = 4.10/Hq which is enhanced over that 
of the ISM, while helium is relatively depleted: y^^'^ = 0.067 = 0.69He/HQ. Thus we use the 
present-day cosmic-ray abundances in the scaling rule for the source particles: 

Xi{t) 



(7) 



This scaling has the desired properties that yi{t) oc Xi{t), and that yp^(t) = yp^(to) when 

Xi{t) = Xi^Q. 



4-1-2. Superbubble Accelerated Particles 

An accelerated particle component in addition to the GCRs has been proposed utilizing i) 
individual massive exploding stars (Casse et al 1995) and ii) collective explosions in superbubbles. 
These are large regions of rarefied, ionized, metal-rich gas associated with star formation regions. 
Observationally, superbubbles are identified by their dense shells of neutral gas. The regions 
interior to these shells are devoid of neutral gas, but sometimes show evidence of ionized gas. 
Superbubbles probably play a fundamental role in the structure and energetics of star-forming 
regions, and of the ISM in general; for a review see, e.g., Spitzer ( p.990| ), Tenorio-Tagle & 
Bodenheimer ( |198^ ), MacLow and Mc Cray (|98|), Walker et al (|l99|). 



The structure and dynamics of superbubbles make them favorable sites for particle 
acceleration. The collective effects of massive stars which sweep material out into the dense shells, 
leaving behind hot, rarefied interiors. The copious injection of matter and energy by massive 
stars provides the ingredients of a powerful ion accelerator. Ejecta from core collapse supernovae 
generate recurrent but weak shock waves, and turbulence. These shocks act as accelerators within 
the superbubble environment. Moreover, these regions provide a large energy reservoir: they have 
the necessary power (lO^'^ ergs and higher) and size to energize a great number of low energy 
particles. 

Superbubble accelerated particles (SAP) generated this way could be a significant a 



Galaxy-wide component of non-thermal particles (Vangioni-Flam et al. 1996, 1998). The SAP's 
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would have a metal-rich composition reflecting that of the superbubbles. According to recent 
calculations by Bykov (1999), the SAP energy spectrum is not very different from that of OCR's 
(but other spectral shapes are possible). The SAP's would inevitably produce BeB over to the long 
duration of massive star activity in the parent OB associations (10^ — 10^ yr), and thus are likely 
to be a significant BeB source. The production of BeB in superbubbles, as well as gamma-rays 
and X-rays (Vangioni-Flam, Casse, & Audouze |1999| ) is thus well- motivated. The existence of 
these particles can be confirmed by gamma-ray line observations that will be performed by the 
European INTEORAL satellite. One could even turn the problem around, and use constraints on 
primary BeB to probe the nonthermal component and physical conditions in superbubbles. 

The study of the particle acceleration in superbubbles is only in its infancy. Much remains to 
be explored, including the relation between the parent molecular clouds and their stellar products, 
the structure and distribution of the environment and its degree of ionization, and the parameters 
of the shocks. All of these will affect the composition and the energy spectrum of the accelerated 
particles. 

The spectrum, (/"sBjint) of energetic particles within the rarefied superbubble interior has been 
studied in considerable detail (Bykov |1999| ). Due to the low densities of these regions, nuclear 
interactions are slow and thus BeB production in the interiors is negligible. The dominant BeB 
synthesis instead occurs as the particles escape from the interiors and traverse the dense shells. 
Thus, in deriving BeB production rates, the interior flux, 0sB,mt) becomes the source for the 
particle flux within the shells surrounding them. That is, QsB.sheii (;/)sB,int- For the (/)sB,int, 
we use a simple analytic form which closely follows Bykov's ( |1999| ) C = 0.1 model. The source 
spectrum appears in Figure Ipi. As one can see, it does not differ substantially from the OCR 
source spectrum. 

The energetic particle composition reflects that of the superbubble interior, and thus is 
highly enriched in supernova products. To model this, we take the particle composition to be 
time-independent in the early Oalaxy. We take the superbubble abundances to be the same as the 
ejecta of massive stars with primordial Z = composition. Specifically, we set the composition to 
be that of a 40Mq star, Woosley & Weaver's ( |1995D model U40B. The particles are thus highly 
enriched (relative to solar) in O {y^ = 320/Hq), and also enriched in He and C. 

4.2. Cosmic Ray Propagation 

We begin with a transport equation which describes the propagation and deceleration of 
nuclei from their sources through the ISM. This equation is a general one, and can be applied 
to both OCR and SAP components. Once launched in space, accelerated nuclei are deviated 
by the magnetic field irregularities and lose memory of their origin. The particles propagate 
diffusively throughout the Oalaxy. As they traverse the ISM, they suffer ionization energy losses 
and undergo nuclear interactions. For the lower energy particles, ionization losses are important 
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and the distances traversed are short. High energy nuclei, however, traverse comphcated paths 
whose length is much larger than the thickness of the Galaxy. Each time they reach the border 
of the Galaxy (badly defined) the particles have a small probability of escape. The Galaxy in 
this context is like a "leaky box," and we will describe cosmic ray propagation via the leaky box 
model. Here we briefly summarize this model, which is discussed in detail in the context of BeB 
nucleosynthesis by, e.g., Prantzos, Casse, & Vangioni- Flam ( |1993| ) and Fields, Olive, & Schr amm 



(|1994|) 



The leaky box propagation equation can be expressed as 

^^N{E) = Q{E) - \n{E) + A N{E)] (8) 

where the last equality assumes a steady state. The source (injection) spectrum Q{E) depend on 
the cosmic ray model, as discussed in the preceding sections. Coulomb interactions with ambient 
electrons lead to an ionization energy loss rate h{E). The escape time, r, is the harmonic mean of 
the confinement time and the nuclear destruction time. The energy distribution is modified in the 
course of the propagation since energy losses are energy dependent. 

Dividing by the mean density of the ISM one gets 

-^ct>{E) = q{E) - + A l^^E) (9) 

with u>{E) = dE/dX, in MeV/(gcm^^). The equilibrium spectrum, 0, is solution of this equation. 
BeB production rates and ratios are obtained by integrating over this spectrum rather than the 
source spectrum q. Propagated spectra for GCR and SAP models appear in Figure l^b. 

The production rate of ^ € BeB per unit volume is 

-ni = J2^iJ dE afj{E) 4>j{E) S,{E) (10) 

Here nj is the ISM number density of species i, 4>j{E) is the solution to eq. (^ for particle species 
j, and crfj{E) is the cross section for the spallation reaction i + j ^ £ + ■ ■ ■ (Read & Viola 1984|) . 



The factor Si^E) = exp[— i?^(£')/A] is the probability that the daughter nucleus £ is stopped and 
thermalized in the ISM before it can escape from the Galaxy. 



4.3. Neutrino Process 

The neutrino process (Woosley et al. 19901) provides a means for direct (primary) ^^B 
production in supernovae. As the large neutrino flux streams from the supernova core, it traverses 
the overlying shells of material. When the flux passes through the carbon shell, some of the 
neutrinos collide inelastically with the ^^C nuclei, and make ^^B spallatively, by removing a 
nucleon: ^^C + v ^ ^^B + p + v' . Other inelastic reactions create ^'^B and ^Be, but these species 
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do not survive the subsequent passing of the blast wave. Thus, the j/-process provides an elegant, 
primary source of ^^B, but does not make Be. 

We adopt the //-process yields of Woosley &: Weaver ( 1995 ), who found that ^^B is made at 
significant levels. The precise yields, however, depend on the neutrino temperature, and hence are 
uncertain. Although the ^^B production is very sensitive to the temperature, the relative yields 
for different mass supernovae are probably better known. To allow for this uncertainty, we follow 
previous work (Olive et al. 1994 ; Vangioni- Flam et al. |1996| ; Fields & Olive |1999aD in taking the 
i/-process yields to be uncertain by an overall factor. We thus write mej,ii(m) = /,ym^^(m), 
where the yield mej,ii of ^^B is a function of progenitor mass m, and scales with the Woosley & 
Weaver yield m^_^(m). The scale factor fy will be set by a fit to the data, as described below, 
and the resulting value amounts to a constraint on the neutrino temperature. 



4.4. Chemical Evolution 

The production of BeB depends on the intensity and composition of both GCR and SAP, as 
well as the ISM composition and gas mass. These quantities are time dependent, and to follow 
their evolution requires a model of Galactic chemical evolution. The chemical evolution model 



adopted here is described in detail in Fields & Olive 1999a, so we will only summarize the key 
points. 

For the star formation rate ilj{t) we adopt = AMgas, with A = 0.3 Gyr^^; our results are 
insensitive to the details of i/j. We use a power law initial mass function, oc with 
m G (O.IMq, IOOMq). For simplicity we use a closed box model, i.e., without infall or outflow. 
We account for different stellar lifetimes as a function of mass, and thus we do not use the 
instantaneous recycling approximation. 

The elements evolved in the code are H, He, C, N, O, and Fe, as well as BeB. The high 
mass (Type II supernova) yields are from Woosley & Weaver ( |1995| ). The intermediate mass 



(AGB) yields are from van den Hoek & Groenewegen ( 1997 ). In both cases, the full metallicity 
dependence of the tabulated yields is used. 

The different roles of O and Fe deserve mention. We stress that the oxygen yields are the 
most critical, since as we have argued O is the better diagnostic for BeB origin, and thus the 
BeB-0 trends are the main focus throughout this paper. The adopted O yields of Woosley & 
Weaver ( 1995| ) are in broad agreement with the work of Thielemann, Nomoto, & Hashimoto 



( 1996 ). Furthermore, the O yields are not strongly altered in explosive nucleosynthesis, nor are 



they sensitive to the mass cut. 

The yields of iron, on the other hand, are much more uncertain. However, because Fe has been 
shown in all previous work, we also will show this for comparison. Our treatment of Fe production 
follows that of Fields & Olive ( 1999a| ) , and is driven by the new O /Fe data. The observed O /Fe-Fe 
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relation (eq. ||) is key to interpreting the BeB-Fe slopes, so we adopt this relation, and solve for 

Fe/Fe© oc (O/O0)^/(^+'^o/F^) (11) 



iron: 



Since O yields and evolution are fixed in the model, this strategy amounts to fixing the effective 
yield for Fe. Note that the non-constancy of 0/Fe (i.e., WQ/Pe 7^ 0) demands that the derived Fe 
yield is not proportional to the O yield. 

In other words, we can derive the effective Fe yields needed to agree with the 0/Fe data and 
the O yields. The effective Fe yield (mej^Pe) is the mean Fe mass ejected per supernova, which is 
the total Fe ejection rate divided by the supernova rate: 

(mej,Fe) = (12) 

iVii 

The iron mass ejection rate can be inferred as follows. The usual chemical evolution 



expressions give Epe = MgXpe + EXpe- We compute Xpe via eq. (|llD , which is imposed via 
Xpe/Xpe = (1 + ujq/pq)~^Xo/Xo, where Xq is computed as usual from the high-mass yields 
within the model. Thus (mej^Pe) is an effective yield, averaged over the mass function, at each 
epoch t. 

Figure |2| shows the time evolution of (mej,Fe), and for comparison, a similar effective yield for 
O. At any give time, the mass range includes all of those stars whose mass m is larger than mt, 
the mass whose lifetime T(m) = t. Thus, the average includes an increasingly larger mass range of 
stars as time goes on. We therefore also plot these results as a function of mt in the lower panel of 
Figure ||. We see that the effective O yield starts high (> IOMq) when the progenitors have high 
mass, and then levels off at IOMq, where the average then includes all supernovae. In contrast, 
Fe starts low and then builds up. Within the m > IOMq regime, the Fe yield is not constant, but 
spans a large range, to include values at lower masses which are of the order of Fe yields estimated 
for observed supernovae, though somewhat lower than the yields of Woosley & Weaver ( 1995| ). 



Above 60 Mq, the Fe yield is assumed to be negligible. Below m < IOMq, (m-ej^Fe) continues to 
rise to values approaching IMq. This rise thus coincides with the epoch when Type la supernovae 
should begin to contribute to iron enrichment. The Type la iron yields are expected to be nearly 
constant per explosion event, and to approach IMq] thus, one can view the m < IOMq behavior 
of (77iej,Fc) as a measure of the increasing ratio of Type la/Type II supernova events. 

It is noteworthy that the inferred Fe yields are at variance with the yields that Shigeyama 
& Tsujimoto ( |1998 ) inferred from Mg/Fe ratio in ultra-metal-poor halo stars (and from the light 



curves of supernovae, which however arise from a different population than the halo stars). While 
the trend in (mej^Fe) is inferred from 0/Fe to be decreasing with increasing progenitor mass, the 
trend inferred from Mg/Fe is found to do just the opposite. This difference is indeed striking 
but not surprising, since Shigeyama & Tsujimoto (|1998|) adopted a Mg/Fe relation that decreases 
towards low metallicity, rather than increases as 0/Fe does. This result discrepancy indicates the 
need to get good observational constraints on the 0-Fe-Mg relations. 
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An alternative model has been proposed in which the 0/Fe relationship is obtained by 
assuming differential mixing of Be, O and Fe, the latter being locked up into grains (Ramaty and 
Lingenfelter 1999, Ramaty et al. [1999 ). With the nominal yields of Woosley and Weaver (1995), 
they obtain an acceptable 0/Fe evolution at the expense of an additional free parameter: the 
mixing time. 



5. Model Results 

In presenting results, we will compare calculated Be and B trends with observations. In 
doing this, we focus on oxygen as a metallicity indicator. Using oxygen rather than iron has 
several theoretical advantages. Oxygen is a good tracer of the cosmic ray accelerators (Type 
II supernovae). Also, oxygen (and to a lesser extent C and N) provides the main raw material 
from which Be and B are carved by spallation. And while both oxygen and iron are both made 
by Type II supernovae, the oxygen yields are considerably less uncertain as they are dominated 
by well-understood hydrostatic burning processes, as opposed to the iron yield which depends 
sensitively on the details of the (poorly understood) explosion mechanism and mass cut. However, 
as noted above, Mg may be a promising evolutionary tracer: from the theoretical point of view, 
Mg shares the same robustness in prediced yields that O does, since it is not affected by the poor 
knowledge of the position of the mass cut between the neutron star and the supernova ejects on 
one hand; from the observational point of view, Mg abundances depends less on the details of the 
stellar atmospheric models (for instance it is less subject to NLTE effects, Fuhrmann ( |l99j )). 



5.1. Be and B Normalization 

The three components in the BeB evolution lead to three free parameters: 

1. The normalization of the GCR component. We scale the GCR flux via 
$(t) = /gcr[V'(*)/V'(*o)]^05 where the star formation rate scaling, 'ip{t)/ip(to), is 
computed self-consistently in the model, and $o is the total flux in present-day cosmic rays, 
as in Figure 1. 

2. The SAP (energetic particle) normalization. We put <I>sAp(i) = fsAP^(t). 

3. The z^-process normalization. The scaling fi, is described in § |4.3| . 

Each of the fi are free, in that there are independent constraints on the normalizations, but we 
do have some independent constraints on their values. For example, the GCR normalization 
is essentially fixed by the present Galactic average GCR flux, which is constrained by GCR 
flux measurements in the vicinity of the earth at high energy (> 1 GeV) as well as by 7-ray 
observations (see Strong and Moskalenko ( 1999[ )). Thus, the GCR normalization is not arbitrary. 
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but on the other hand, since the present-day Galactic average flux is uncertain to within a factor 
of a few (Mori 1997). We find that the scaling is /gcr = 2.1, well within the uncertainties of the 
Galactic average GCR flux. Note that the other two normalizations are more poorly constrained; 
indeed, the BeB data are probably the best way to constrain them. 

We thus need three pieces of data to fix the normalizations; then the resulting Be-B-0 trends 
are predictions. This is done as follows 



1. A break point [0/H]eq is chosen to match that of the data set being fit. This fixes the 
relative strength /sap of the GCR vs SAP components. Thus, the evolution of the species 
that have only GCR and SAP sources, namely Be and ^'^B, are fixed up to an overall scale 
factor. 

2. The Be evolution is normalized to reach (Be/H)© at (0/H)q, which fixes /gcr- Thus Be-0 
is forced to have the right [0/H]eq, and must go through the solar point.^ We have now fixed 
the GCR and SAP components completely, only using the Be break point and solar point. 
In doing this, ^"B is completely determined, and ^^B is determined up to the z^-process 
contribution; both of these are predictions of the model. 

3. Finally, the i^-process contribution to ^^B is scaled so that the predicted ^^B/^'^B ratio agrees 
with the meteoritic value at solar metallicity; this fixes f^- Now the model is completely 
determined, and the scalings with O cannot be adjusted further. 



Note that this procedure does not force elemental B to go through the solar point. As we will 
see, the models nevertheless do pass through or very near the solar point, which should be viewed 
as a success of the scenarios being tested. 



5.2. Comparing Theory and Observation 

To compare theory and observation we have run models following the procedure just described, 
and we now plot them with the observations for the different sets of stellar parameters. 

As discussed above, the break point [0/H]eq is an input to the model. However, it is not an 
entirely free parameter, as is evident when applying a two-component (linear plus quadratic) fit 
to the Be-0 points from the GCR model alone. We find that this model has a small but nonzero 
linear component, ai = 0.043, in addition to the expected large quadratic term, 02 = 3.8 value. 
The small linear component results from the fact, shown in Fields & Olive ( |1999a ), that the Be-0 



slope is slightly less than 2, due to small astration effects. This flattening leads to a small but 



^ Fortunately, the solar abundances of Be and B are now well established: there is now good agreement between 
the photospheric and meteoritic abundances of both beryllium and boron. After a long debate, the situation has 
finally been settled: Be and B are undepleted in the solar photosphere. 
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nonzero linear component in the fit. Because ai > 0, the GCR model alone has a break point at 
— 1.94, without any addition of the primary process. Thus, there is no way to arrange a break 
point below —1.94, since adding a purely primary component only raises the value of ai and thus 
the break point, defined by 01/02- 

In making the comparison between theory and observations, we have chosen the break point 
for each method corresponding to the lowest number of data points. This data subset uses stellar 
atmospheric parameters which have been derived directly from the prescription discussed, and 
therefore should be most free of systematic effects. For example, in the Balmer case, only 19 of 
the total 31 stars are truly what we have called Balmer. Recall an additional 8 were described 
as Balmer + King in section 2 and may have additional systematic effects entering which we are 
trying to avoid. This intermediate set (as well as the full data set) were used to compute slope 
values, etc., for comparison only. 

The Balmer data have [0/H]eq = —1.75 as given in (Table |5|); the results of a BeB model with 
this [0/H]eq value appears in Figure ^. As it happens, this break point is close to that of the GCR 
component, and thus there is only a small SAP component for this model. Note that in this figure 
and those that follow, the SAP component and the GCR components cross at value of [0/H] below 
[0/H]eq. This is due to the small primary component in the GCRs themselves (oi^gcr = 0.043), as 
discussed above. Thus [0/H]eq is determined by the sum of the two primary components, relative 
to the secondary component. One can see that the Be-0 trend in this model fits the data quite 
well. The B-0 model contains in addition the z^-process. The lower panel of Figure ^. plots the 
different components as well as the evolution of total abundance which sums the components. We 
see that B-0 is dominated by the primary z^-process at low metallicity, and by the secondary GCR 
production at high metallicity, as expected. However, even at present the z^-process component is 
not negligible, being lower than the total GCR component by about a factor of 2 (as is indeed 
demanded by the fit to the ^^B/^^B ratio). We have scaled the i^-processyields of Woosley et 
al. ( |1990| ) downward by a factor of fi, = 0.4 here and, as it turns out, in the other models as well. 
This scaling also keeps the production of ^Li small enough so as not to affect the Spite plateau as 



shown in Vangioni-Flam et al. ( 1996 ). In Figure 3, and in the figures that follow, the point shapes 
indicate which data sample the stars fall into (e.g., pure Balmer, Balmer + King, or full). 

Figure Q plots Be and B versus O for IRFMl, which has a break point [0/H]eq = —1.79. 
The upper panel shows Be-0, and is very similar to the Balmer case discussed above due to the 
similarity in the break points. Again, the intrinsic linear portion of the GCR component influences 
the break point, which is why the break point is higher than the metallicity at which the GCR and 
SAP components are equal. The quality of the Be-0 fit is clearly good. In the case of B-0 (lower 
panel), the i/-process contributes as well, leading to a GCR dominance at a higher metallicity than 
for Be-0. Unlike Be-0, the B-0 fit now seems to be a poor one. The overall shape is fine, but the 
curve is clearly high. 

Figure ^ plots BeB-0 evolution for IRFM2, which has [0/H]eq = —1.37. In this case, the 
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transition from SAP dominance at low metallicity to GCR dominance at high metalhcity is 
evident. The results here are qualitatively similar to the IRFMl case: the Be-0 evolution is 
satisfactory, while the B-0 evolution seems high. 

In the next set of figures, we show plots of B/Be versus O. At the lowest metallicity in the 
Balmer case. Be is secondary while B is dominated by the i/-process and is thus primary. Hence, 
the B/Be ratio changes as B/Be oc and has a log slope of -1, as seen in Figure ^. Interestingly, 
even with this rise, the model actually falls short of the two highest B/Be ratios, though the 
qualitative behavior is correct. At high metallicity, the GCR component dominates both Be and 
B, and the ratio flattens. In the IRFM case, at very low metallicity, both Be and B are dominated 
by the SAP component, and again the curve flattens. The behavior of B/Be for the other two 
cases are shown in Figure (for IRFMl) and in Figure | (for IRFM2). In the latter case the 
break point is higher and thus there is only a narrow range in [0/H] where B/Be drops, and so 
the overall change in B/Be is smaller. 

In the IRFMl case, B/Be is also far from the very high B/Be data points at the lowest 
metallicities. This illustrates the power of the B/Be ratio to constrain the models. The data is 
still poor however. While the situation is presently too uncertain to clearly exclude the models 
presented, it is clear that more B/Be data at low metallicity will play a decisive role in determining 
Be and B origins. 

In Table 0, we show the production ratios of GCR and SAP at solar metallicity. That these 
are so close traces back to the very similar energy spectra (in spite of their different composition). 

Table 7: Production ratios for GCR and superbubble accelerated particles (SAP) 



Ratio 


Present (Z = Zq) GCR 


SAP 


^Li/^Be 


5 


5.8 


iOB/9Be 


5 


4.9 


7Li/6Li 


1.4 


1.2 


^Li/^Be 


7 


6.8 


iiB/9Be 


12 


11.7 


iiB/iOg 


2.5 


2.4 


Li/Be 


12 


14 


B/Be 


17 


16.6 



5.3. Energetics 



As emphasized by Ramaty et al. ( 1997 , 199S| ), a study of the energy budget associated with 
BeB production models can give additional insight. There are several ways to proceed, but the 
key point is to link the energy in accelerated particles required to make the BeB with the available 
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energy budget in the form of supernovae. The power of the analysis is in judging whether the 
energy requirement changes with metalhcity (and thus epoch); in particular, it is important to see 
whether there is too large an energy requirement in Population II. 

We find that there is not an energetics problem in the models we consider. One way to see 
this is as follows. Given a model, we wish to compute the accelerated particle energy input per 
supernova: 

^CR 

eCR = (13) 

^VsN 

where <?cr(*) is the total energy going into particle acceleration throughout the Galaxy per 
second, and A^sN is the supernova rate. The direct approach is simply to note that in both 
the GCR and primary components, supernovae are the power source of energetic particles: 
"^CR oc Q oc A^'sN. Since the power input is proportional to the supernova rate, their ratio ecR 
is constant, independent of metallicity. The models always require the same accelerated particle 
energy per supernova. Note that this approach simply examines the internal consistency of the 
model. We have not used the properties of the Be data. 

We can use the data, however, and get the same conclusion. One can derive the ecR scaling 
observationally by relating the Be-0 and Be-Fe data to the energy budget. This comes about by 
using the definition of ecR above (eq. |l^) and a "chain rule" expansion: 

VMBc/inodcl V^^/obs V^SN/j„odcl 

for the metal tracers Z G 0,Fe. As we have indicated, the first and last terms must be taken from 
the model (note that in Pop II, Scr/Mbb — W/Q{Be) in the Ramaty et al. 1997| notation). The 



middle term can be related to the data, since M-Q^jMz oc Be/Z as long as Be can be expressed 
as a power law in Z . All theory terms (both denominators and numerators) are individually 
calculable, so we can now rearrange terms: 

V^SN/^odel V^^^/obs V^^Be/jnodei 
(Be/Z)obs , . 

°^ m I7\ (1^) 

(Ce/Z jmodel 



Equation (16) follows from eq. (15) due to the intrinsic constancy of "fcRZ-^SN in the models; 
now, however, something new has been added in that the observed Be/Z ratio appears explicitly. 
In particular, we see that the scaling of ecR with metallicity depends on the ratio of observed to 
predicted Be/Z. If the observations are well matched by the predictions, then the ratio remains 
constant, as does ecR- In other words, the constancy of ecR is determined by the goodness of the 
Be-Z fit. Thus, all models which fit the data are able to avoid an energetics problem. This applies 



quite generally, to the models presented here as well as others (e.g., Ramaty et al. |1999D which fit 
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the data.0 

This argument shows that the models in question are energeticahy sound, simply due to their 
agreement with the data; however, a more detailed analysis is instructive. As we have argued, it 
is preferable to use the O data as the metal tracer Z, as this is the most direct diagnostic of BeB 
origins. As we have emphasized, and shown in Table ||, the observational data clearly show that 
Be/0 is not constant, but rather (Be/O)obs ~ o*^-^^*^'^, where the exact power depends on the 
atmospheric model used. This factor alone would lead to a decrease in ecR at low metallicity. 
However, when the GCR is dominant over SAP's, (i.e., for metallicities [0/H] > [0/H]eq), the 
"production efficiency" factor scales as <^cR/-^Be ~ A''sn/(OA''sn) ~ O^"*^. This reflects the reduced 
efficiency for GCR Be production in the early Galaxy: the paucity of O targets in the early ISM 
leads to a lower BeB production for each supernova's complement of accelerated particles; this 



factor alone would lead to an increase in £cr- Thus, the first two factors in eq. (14) offset each 
other, which is largely why ecR is constant. 



We must consider several smaller and more subtle effects to complete our analysis of eq. (14). 
The last factor in eq. (0), the "average yield" Mz /^SN has a weak dependence on O, rising 
slightly towards low metallicity. However, we have until now assumed that M-q^/Mz oc (Be/Z)obs- 
This is a good approximation in Pop II, where Be evolution is dominated by production only, i.e., 
-^Be ~ Q(Be). But in Pop I, astration effects become important since many low-mass stars begin 
to die and return their astrated. Be-free gas, which dilutes the ISM Be abundance. Astration 
manifests itself in the Be-0 and Be-Fe relations by the turnover of Be at high metallicity: Be 
production is overtaken by destruction. Thus, for Pop I, the change in Be gas mass is lower than 
the cosmic ray production term would give: Mbc < (5(Be). Pop I Be production is, in this respect, 
less efficient than in Pop II where astration is small and Mbc — Q(Be). Thus, the Be-Z relation 
is flattened in Pop I, and the extrapolation of ecR via eq. ( p!4| ) back to Pop II overestimates the 
energy needed if this astration is neglected. Finally, at very low metallicity, [0/H] < [0/H]eq, the 
SAP component takes over the energy budget. Thus, just as SAP's are primarily responsible for 
BeB production below [0/H]eq, they are also the dominant energy source in this range. 



Finally, we consider using Z = Fe in eq. (|T4|). To reiterate, the question regarding the 
energetic problem, rests on ecR- This quantity may be constant, even though W/QiJie) is not 
constant, depending on the observed Be/Fe ratio and the ratio of Be/Fe produced. Since the 
Be-Fe slopes are larger than (though close to) 1, Be/Fe increases with Fe. More significantly, the 
Be/Fe production ratio does not scale as Fe itself, but rather Mbg oc O and O ~ Fe*^'^. Similarly, 
the iron yield per supernova may also not be constant if the yield of Fe falls off at large stellar 
masses (see Fig ^). Both of these facts simply reflect the fact that 0/Fe is not constant in Pop II. 
Thus, models in which W/Q{Be) is not constant, can not a priori be disregarded. 



^Strictly speaking, the scaling argument we have given says that models which do not have an energetics problem 
at the present epoch and fit the Be/Z data will also satisfy the energy requirements in earlier epochs. 
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6. Conclusions 

The history of the accelerated particles within the Galaxy is coded in the evolution of 
spallogenic Be and B. Recent theoretical studies of particle acceleration suggest that the Galaxy 
may have two components of nonthermal particles. These two populations yield different BeB 
evolution, GCR particles lead to secondary (oc O^) production, while nuclei accelerated in 
superbubbles give a primary (oc O) production. 

We have searched for the corresponding evidence of two components in the evolution of Be 
and B versus O and Fe abundances in low metallicity halo stars. Specifically, we have made 
two-component fits of abundance data, to identify the presence of primary and secondary BeB 
components, and to quantify their relative strength. This analysis has been carried out using 
consistent stellar atmosphere parameters based on Balmer lines and the infrared flux method 
(IRFM). Some results depend on the method used to get the atmospheric parameters, while others 
do not. In both Balmer and IRFM methods, the secondary BeB component is seen to dominate 
for [0/H] > —1.4. However, the status of the primary component depends on the atmospheric 
parameters chosen. When the Balmer treatment is used, a primary mechanism is not required 
over the observed metallicities, whereas it is required in the IRFM2 case, where it dominates at 
[0/H] < -1.4. 

We use the BeB-OFe fits to constrain the theories, whose predictions we model by combining 
cosmic ray nucleosynthesis and Galactic chemical evolution. The GCR component has the 
observed present-day spectrum, is assumed to be accelerated from supernovae, and to have a 
composition which scales with that of the ISM. For the case of superbubbles, the energy spectrum 
has been taken from Bykov (1999), and is quite similar to the that of the GCRs. Thus there a 
little hope to differentiate the two components on an isotopic basis, since the production ratios of 
superbubbles and GCR are very similar in spite of their different composition (C/0 = 0.09 against 
C/0 = 0.81 respectively). The two-component BeB-0 fits constrain the models by demanding 
that: (1) the GCR production dominates over the superbubble component at [0/H] > —1.4 in 
all cases, and in the Balmer case it is unclear whether the SAP component is needed; (2) the 
secondary, GCR component requires that the mean GCR flux at present factor of 2.1 higher than 
current estimates; given the uncertainties, this is quite plausible; and (3) the i^-process is needed 
to produce ^^B, but at a level that is 40% of the flducial Woosley & Weaver ( |1995 ) yields; again, 
this is well within the uncertainties of the yield calculations. 

Our results can be further tested, and the constraints tightened, by several types of 
observations. Simultaneous and high signal/noise measurements of Be, B, O, and Fe, are called 
for. These must be interpreted within detailed and consistent atmospheric models including 
accurate NLTE corrections. The current situation is not definitive, since the conclusions regarding 
the primary component depend sensitively on the atmospheric parameters adopted. A critical goal 
is to determine Be and B abundances-as well as key ratios such as B/Be, Li/Be, and ^^B/^'^B-at 
very low metallicity, where the primary and secondary origins differ most in their predictions. 
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It is just as important to firmly establish the nature and evolution of 0/Fe; other metallicity 
indices, such as Mg, could help to clarify this question. These should be analyzed systematically 
as we have tried to do with oxygen. Gamma-ray line astronomy should also help to confirm the 
possibility of acceleration of freshly synthesized nuclei in superbubbles and OB associations. This 
will be a key objective for the European INTEGRAL mission. Thus, we are pleased to report that 
the nature of LiBeB evolution, is a problem that is ripe for new and precision observations, which 
can go far to reveal the origin and history of accelerated particles in the Galaxy. 

Wc thank Andrei Bykov for graciously making available data files for his published spectra. 
We thank Roger Cayrcl and Evan Skillman for fruitful discussions. This work was supported in 
part by DoE grant DE-FG02-94ER-40823 at the University of Minnesota. 
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Figure Captions 

1. Source and propagated spectra for accelerated particles. 

Top Panel: Source spectra Q{E) for the GCR (solid curve) and SAP (dashed curve) 
components. E is the kinetic energy per nucleon in MeV/n, while the units of Q{E) are 
arbitrary, and the relative scaling of the two components has been adjusted for clarity. 
Bottom Panel: Propagated spectra corresponding to the sources in the top panel. 

2. Oxygen (solid line) and iron (dashed line) yield estimates, calculated from the numerical 
model. For derivation see discussion in text. 

Top Panel: Yields versus time. 

Bottom Panel: Results from top panel translated to a mass scale via the stellar mass-lifetime 
relation T(ni). 

3. Be vs O {top panel) and B O {bottom panel). Data shown are the Balmer points, which 
are found (Table ^) to have a break point as indicated. Models are adjusted to have the 
break point and O /Fe slope of these data. 

4. As in Fig. |3|, for the IRFMl data and corresponding model. 

5. As in Fig. |3|, for the IRFM2 data and corresponding model. 

6. The B/Be ratio vs O. Shown are the Balmer data and the corresponding model. 

7. As in Fig. |6|, for the IRFMl data and corresponding models. 

8. As in Fig. |6|, for the IRFM2 data and corresponding models. 

9. Be and B vs Fe, for the Balmer data and break point. As discussed in the text, Fe in the 
model is obtained by scaling from O, using the observed relation [0/H] = WQ/PePe/H], with 
coq from Table |l|. 

10. As in Fig. |9|, for the IRFMl data and corresponding models. 

11. As in Fig. |9|, for the IRFM2 data and corresponding models. 
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